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1  Schematic  diagram  of  a  2-Dimensional  Wolfhard  -  Parker  diffusion  flame  burner, 

showing  burner  geometry.  The  burner  has  been  modified  to  introduce  mist  into  the 
air  co-flow  through  mist  generation  chambers  located  next  to  the  oxidizer  channel. 
A  fraction  of  the  mist  generated  in  this  chamber  is  entrained  into  the  air  stream  and 
interacts  with  the  diffusion  flame . 

2  Comparison  of  numerical  and  experimental  results  for  base  case  CHa  -  Air  diffusion 

flame.  Numerically  computed  temperature  profiles  are  compared  with  thermocouple 
temperature  measurements  at  a  height  of  5,  15,  30  and  55  mm  above  the  burner 
surface . 

3  Comparison  of  numerical  and  experimental  results  for  CH4  and  18.2%02,81.8%N2 

diffusion  flame.  Numerically  computed  temperature  profiles  are  compared  with  ther¬ 
mocouple  temperature  measurements  at  a  height  of  6,  40,  46  and  50  mm  above  the 
burner  surface . 

4  Comparison  of  numerical  and  experimental  results  for  CH4  and  16.5%02,83.5%iV2 

diffusion  flame.  Numerically  computed  temperature  profiles  are  compared  with  ther¬ 
mocouple  temperature  measurements  at  a  height  of  6,  26,  46  and  60  mm  above  the 
burner  surface . 

5  Contours  of  temperature  and  heat  release  rate  over  a  CH4  —  Air  diffusion  flame 

burner.  The  fuel  jet  velocity  is  2.81cm/sec  and  air  jet  velocity  is  25.87cm/sec.  The 
flame  height  is  approximately  7cm  and  the  maximum  flame  temperature  is  1970A'. 
The  heat  release  rate  profile  shows  the  region  of  exothermic  chemical  activity  above 
the  burner  surface . 

6  Contours  of  methane  CH4  and  oxygen  O2  density  above  a  CH4  -  Air  diffusion 

flame  burner.  The  fuel  jet  velocity  is  2.81cm/sec  and  air  jet  velocity  is  25.87cm/sec. 
These  calculations  are  performed  using  a  single  step  arrhenius  chemistry  to  study 
methane-air  combustion . 

7  Contours  of  product  species  carbon-di-oxide  CO2  cind  water  vapor  H2O  density 

above  &CH4-  Air  diffusion  flame  burner.  The  fuel  jet  velocity  is  2.81cm/sec  and 
air  jet  velocity  is  25.S7cmf  sec.  These  calculations  are  performed  using  a  single  step 
arrhenius  chemistry  to  describe  methane-air  combustion . 

8  Velocity  vectors  and  stream  line  pattern  during  methane-air  combustion  above  a 

Wolfhard  Parker  diffusion  flame  burner.  The  velocity  vectors  have  been  color  coded 
with  red  showing  the  highest  gas  speed  and  blue  showing  the  lowest  gas  velocity.  The 
streamline  pattern  (black)  clearly  indicate  the  entrainment  of  oxidizer  along  with  the 
suppressants  into  the  diffusion  flame.  The  fuel  jet  velocity  is  2.81cm/sec  and  air  jet 
velocity  is  25.87cm/sec . 

9  Suppression  of  diffusion  flames  subjected  to  50;u  diameter  fine  water  mist  droplets 

injected  with  a  velocity  of  25cmf sec  into  the  air  flow.  Figure  shows  reduction  in 
flame  temperature  and  heat  release  rate  profiles.  The  temperature  and  heat  release 
rate  profiles  indicate  an  increase  in  flame  height  and  flame  suppression  due  to  sub¬ 
extinguishing  application  rates  of  water  mist.  The  heat  relecise  rate  profile  also 
indicates  the  presence  of  an  endothermic  heat  release  profile  (-ve  values)  showing 
the  approximate  location  of  evaporation  of  the  droplets . 
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10  Suppression  of  diffusion  flames  subjected  to  150p  diameter  fine  water  mist  droplets 

injected  with  a  velocity  of  250cm/sec  into  the  air  flow.  Figure  shows  reductions  in 
flame  temperature  and  changes  to  the  heat  release  rate  profile  as  compared  to  the 
base  case.  The  temperature  and  heat  release  rate  profiles  indicate  an  increase  in 
flame  height  and  flame  suppression  due  to  the  sub-extinguishing  application  rates  of 
water  mist.  The  heat  release  rate  profile  also  indicates  the  presence  of  an  endothermic 
heat  release  profile  (-ve  values)  showing  the  approximate  location  of  evaporation  of 
the  droplets . 

11  Suppression  of  diffusion  flames  subjected  to  50/i  water  mist  droplets.  Figure  shows 

temperature  contours  for  various  spray  injection  densities  of  1000,2000,3000  and 
4000drops/cm^.  These  droplets  are  injected  along  with  the  air  with  an  initial  injec¬ 
tion  velocity  of  25cm/sec.  Also  shown  are  the  bcise  case  profiles  indicating  various 
levels  of  suppression . 
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10  -  20//,  20  -  30//,  30  -  40//  and  40  -  50//.  Droplet  section  densities  for  three  of 
the  five  sections  have  been  shown  to  illustrate  the  location  of  maximum  density  and 
position  where  these  droplets  evaporate  to  form  smaller  and  smaller  droplets  (sections). 
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of  the  five  sections  have  been  shown  to  illustrate  the  location  of  maximum  density 
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19  Temperature  contours  above  a  methane-air  diffusion  flame,  a)  Base  case  contours, 

b)  Temperature  contours  during  the  injection  of  50/i  water  droplets  with  a  spray 
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Numerical  Modeling  of  Fire  Suppression  using  Water  Mist. 
1.  Gaseous  Methane- Air  Diffusion  Flames 


1.  INTRODUCTION 

Water  as  a  means  of  fire  suppression  has  been  in  use  from  ancient  times.  The  phase-out  of 
halons  and  search  for  alternative  technologies  that  preserve  all  of  the  benefits  of  a  clean  total 
flooding  agent  without  the  adverse  environmental  impact  has  sparked  renewed  interest  in  water 
mist  technology.  The  need  for  low  weight  impact  replacement  sprinkler  systems  on  commercial 
ships  has  been  driven  by  International  Maritime  Organization  regulations  requiring  retrofit  of  fire 
suppression  systems  on  most  commercial  marine  vessels.  This  has  also  given  recent  impetus  to  the 
development  of  low  water  demand,  high  efficiency  mist  systems  to  replace  sprinkler  systems. 

Fine  water  mist  relies  on  relatively  small  (less  than  200  fj.m)  droplet  sprays  to  extinguish  fires. 
The  mechanisms  of  extinguishment  include  gas  phase  cooling  (thermal  effect),  oxygen  displacement 
by  steam,  wetting  of  fuel  surface  and  attenuation  of  radiative  heat  transfer.  Although  the  usefulness 
of  water  mist  fire  suppression  systems  has  been  demonstrated  in  a  wide  range  of  applications  and  by 
numerous  experimental  programs  (References  [1]-  [14]),  a  widely  accepted  critical  concentration  of 
water  droplets  required  to  extinguish  a  fire  is  yet  to  be  determined.  Factors  that  contribute  to  the 
success  or  failure  of  a  water  mist  system  for  a  particular  application  include  droplet  size,  velocity, 
spray  pattern  geometry,  momentum  and  mixing  characteristics  of  the  spray  jet,  geometry  and  other 
characteristics  of  the  protected  area.  At  this  time  the  effect  of  these  factors  on  system  effectiveness 
is  not  well  known.  There  is  no  current  theoretical  basis  for  the  selection  of  spray  characteristics 
and  other  important  water  mist  system  parameters. 

1.1  Literature  Review 

Extinguishment  of  jet  diffusion  flames  by  the  introduction  of  water  sprays  at  the  base  of  the 
flame  has  been  studied  by  a  number  of  investigators  [Ij-  [5].  These  flames  can  be  extinguished 
very  efficiently  by  this  technique  since  the  normal  air  entrainment  process  assures  that  all  the 
water  spray  added  adjacent  to  the  flame  is  actually  transported  into  the  flame.  The  water  required 
for  extinguishment  is  reported  in  terms  of  the  ratio  of  the  water  supply  rate  to  the  fuel  supply 
rate  at  extinction.  These  ratios  range  from  1.5  to  10.  While  no  systematic  evaluation  of  droplet 
sizes  on  the  water/fuel  ratio  at  extinction  has  been  performed,  the  available  data  indicate  that 
the  ratio  is  reduced  with  decreasing  droplet  size  for  laboratory  flames.  Large  scale  flames  could 
be  extinguished  with  a  water/fuel  ratio  of  1.6.  However,  ratios  of  up  to  10  may  be  required 
depending  on  the  orientation  and  geometry  of  the  spray  nozzles.  While  drop  size  information  was 
not  given,  the  nozzles  used  in  the  experiments  were  expected  to  produce  droplets  on  the  order 
of  1000/im  diameter.  Similar  results  were  obtained  for  laboratory  flames  with  droplet  sizes  of 
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15/iTO.  These  studies  also  indicate  the  reduction  in  flame  temperature  which  occurs  in  flames  with 
sub-extinguishing  application  rates  of  water. 

The  available  research  in  water  spray  fire  extinguishment  was  reviewed  by  Rashbash  [6]  and 
by  Tatem  et  al.  [7].  Rashbash  et  al.  [8]-  [9]  described  two  basic  mechanisms  for  extinguishment 
of  liquid  fuel  fires  :  gas  phase  flame  extinction  and  surface  cooling.  Gas  phase  extinction  was 
observed  to  be  very  rapid  while  surface  cooling  was  much  slower  due  to  the  need  for  water  to  reach 
the  liquid  surface  and  cool  the  liquid  to  below  its  firepoint  temperature.  Fundamental  aspects  of 
solid  fuel  fire  extinction  by  water  have  been  addressed  by  several  researchers  [10]-  [12].  These  works 
have  focussed  on  the  effect  of  water  application  to  the  fuel  surface  and  the  need  for  development 
of  efficient  water  based  systems  for  suppressing  fires.  McCaffrey  [3]  has  reviewed  the  available 
data  and  appropriate  literature  concerned  with  the  application  of  water  sprays  as  a  jet  diffusion 
flame  suppression/extinguishment  agent.  Small  pneumatic  atomizing  nozzles  using  H2  gas,  both 
as  the  flame  source  as  well  as  the  atomizing  driver,  have  been  used  to  scale  high  momentum  jet 
flames  and  to  study  the  effect  of  water  on  the  flame.  The  effect  of  flame  temperature  reduction 
due  to  water  sprays  has  been  observed  to  correlate  with  a  single  spray  parameter-the  median  drop 
diameter.  The  effect  of  adding  water  spray  to  the  gas  flow  below  the  base  of  a  lifted  flame  is  to  shift 
or  raise  the  flame  above  its  normal  position  and  to  lower  peak  flame  temperature  and  radiation 
levels.  Extinguishment  near  blow-off  was  observed  to  be  due  to  shift  in  flame  position.  Though 
McCaffrey  did  not  perform  such  systematic  measurements  for  methane  flames,  the  general  trends 
are  expected  to  be  similar  to  those  of  hydrogen  flames.  It  was  observed  that  for  hydrogen  flames 
that  the  maximum  flame  temperature  per  unit  mass  of  water  per  unit  mass  of  fuel  was  reduced 
from  150  to  50G  as  the  water/fuel  ratio  increased  from  0  to  10.  The  degree  of  cooling  was  also  a 
direct  function  of  the  drop  size  with  higher  cooling  effects  associated  with  smaller  droplets. 

There  has  been  very  little  theoretical  or  experimental  work  to  develop  an  understanding  of  gas 
phase  extinguishment  of  fires.  The  earliest  study  of  the  interaction  of  flames  and  water  mist  was 
performed  by  Seshadri  [15].  He  studied  counterflow  heptane  and  methanol  flames  with  an  oxidizer 
stream  made  up  of  oxygen,  nitrogen  and  water  mist.  This  work  was  of  a  preliminary  nature  in 
that  only  a  single  water  mist  flow  rate  and  drop  size  distribution  was  used.  Based  on  extinction 
oxygen  mass  fractions  as  a  function  of  strain  rate,  he  deduced  one  step  kinetic  parameters  for 
these  systems.  Based  on  the  similarity  of  the  kinetic  parameters,  he  concluded  that  water  had 
only  a  thermal  effect.  No  detailed  flame  measurements  or  droplet  sizes  were  reported.  Lentati  and 
Chelliah  [16]  have  studied  the  effect  of  sodium  bicarbonate  particle  size  in  extinguishing  heptane- 
air  flames.  The  dynamics  of  highly  sooting  fires  in  unbounded  domains  have  been  investigated  by 
Mell  et  al.  [17].  McGratten  and  Stroup  [18]  have  numerically  investigated  conditions  under  which 
vents  and  draft  curtains  are  beneficial,  and  under  which  they  are  detrimental,  to  the  performance 
of  a  sprinkler  system  in  large  enclosures.  Parametric  studies  [19]  have  been  performed  to  predict 
the  actual  delivered  densities  of  early  suppression  fast  response  sprinklers  in  heptane  spray  fire 
scenarios. 

Ewing  et  al.  [21]-  [22]  have  evaluated  the  role  of  thermal  processes  in  the  effectiveness  of  a  wide 
range  of  extinguishing  agents.  Their  work  included  experimental  investigations  of  dry  chemical 
agents  with  critical  particle  sizes  ranging  from  20  -  50//m.  They  observed  that  particle  sizes  below 
the  critical  size  were  fully  decomposed  in  the  flame  while  larger  particles  were  less  effective  since  they 
were  not  fully  decomposed  in  the  flame.  The  work  of  Ewing  clearly  shows  that  dry  chemicals  are 
on  a  weight  bases  the  most  effective  fire  extinguishing  agents  available.  They  found  that  condensed 
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phase  agents  with  large  heats  of  gasification  (evaporation  or  decomposition)  would  be  expected 
to  be  very  effective  agents.  They  proposed  that  heat  extraction  from  the  combustion  region  of 
the  flame  sheet  is  critical  to  flame  extinction.  They  divide  the  flame  into  a  preheat  region  and  a 
combustion  region  and  have  shown  that  flame  extinction  is  governed  by  heat  extracted  from  the 
combustion  region.  Since  dry  chemical  agents  and  water  droplets  can  actually  penetrate  the  preheat 
region  and  are  largely  decomposed  or  evaporated  in  the  combustion  region,  these  agents  are  more 
effective.  While  these  results  provide  interesting  and  potentially  fruitful  directions  for  research  and 
development,  there  is  a  definite  need  for  additional  detailed  experimental  and  computational  work 
to  develop  an  understanding  of  the  interactions  of  water  mist  droplets  with  a  flame  and  the  detailed 
processes  involved  in  extinction. 

Numerical  studies  with  water  mist  in  a  premixed  stagnation  point  methane-air  flame  have  been 
reported  by  Chen  and  Rogg  [23].  Chen  et  al.  [24]  also  presented  numerical  results  for  a  counterflow 
methane/air  diffusion  flame  with  heptane  droplets  added  to  the  fuel  stream.  As  expected  the 
volume  fraction  of  water  mist  at  extinction  is  a  function  of  the  strain  rate.  The  temperature  in 
the  flame  was  a  direct  function  of  the  mist  concentration  and  the  droplet  evaporation  time  was 
predicted.  Atreya  [25]  has  observed  similar  results  for  methane  counterfiow  diffusion  flames.  In 
addition,  he  found  that  the  burning  rate  is  enhanced  at  low  water  addition  rates  and  is  ultimately 
decreased  and  the  flame  extinguished  at  higher  application  rates. 

The  overall  objective  of  this  study  is  to  develop  and  apply  a  numerical  model  for  studying  the 
underlying  processes  involved  in  suppression  of  fires  with  water  mist.  This  report  presents  numerical 
results  on  the  interaction  of  water  mist  with  diffusion  flames  along  with  experimental  validation 
of  the  numerical  model.  The  model  is  then  used  to  understand  the  relative  contribution  of  the 
various  mechanisms  of  extinction  of  flames  subjected  to  a  fine  water  mist.  The  model  provides  a 
detailed  understanding  of  droplet  dynamics  in  a  2-D  flow  field  for  the  study  of  flame  suppression, 
extinguishment  and  mist  entrainment  into  a  diffusion  flame.  The  degree  of  mist  entrainment  into  the 
flame  and  its  dependence  on  droplet  diameter  and  spray  injection  velocity  is  determined.  Numerical 
simulations  have  been  performed  to  evaluate  the  decrease  in  flame  temperature  as  a  function  of 
droplet  number  density.  Droplet  trajectories  are  used  to  identify  the  regions  of  the  flame  (preheat 
zone,  combustion  zone  or  plume  zone)  where  the  droplets  evaporate  and  absorb  energy.  The  droplet 
diameter  and  the  velocity  of  the  droplet  required  to  penetrate  the  flame  under  a  local  application 
scenario  are  also  determined.  Most  importantly  the  model  is  used  to  determine  the  amount  of 
water  required  for  extinguishment,  reported  in  terms  of  the  ratio  of  the  water  supply  rate  to  the 
fuel  supply  rate  at  extinction,  as  a  function  of  droplet  diameter.  Both  local  application  scenarios 
and  a  quasi  total  flooding  scenario  have  been  considered. 

The  next  section  describes  the  governing  equations  that  have  been  used  to  study  the  flow 
above  a  2-Dimensional  Wolfhard  Parker  burner  and  the  interaction  of  water  mist  with  the  diffusion 
flame.  We  have  considered  a  two-continuum  formulation,  wherein  the  gas  properties  and  the  droplet 
properties  are  described  by  equations  in  the  Eulerian  form.  The  next  section  also  describes  the 
single  step  reaction  model  that  have  been  used  to  simulate  the  chemical  reactions  between  the 
fuel  and  oxidizer  species  and  a  sectional  approach  that  has  been  used  to  predict  the  movement 
of  the  fine  water  droplets  along  with  The  numerical  schemes  used  for  time  marching  solutions 
of  the  governing  equations  (gas  phase  and  the  liquid  phase)  along  with  the  initial  and  boundary 
conditions  used,  have  been  briefly  discussed.  The  section  dealing  with  “Results  and  Discussion” 
describes  the  various  simulations  and  parametric  variations  that  were  performed  for  studying  the 
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effect  of  water  mist  on  diffusion  flames.  The  focus  of  this  section  is  on  the  effect  of  droplet  size  on 
suppression  of  diffusion  flames  and  the  relative  contribution  of  the  various  suppression  mechanisms. 
The  experimental  setup  is  discussed  and  validation  of  the  numerical  model  with  the  experimental 
results  has  also  been  provided.  The  degree  of  mist  entrainment  into  the  flame  and  its  dependence 
on  droplet  diameter  and  spray  injection  velocity  is  determined.  We  conclude  by  describing  the  key 
accomplishments  of  this  work  along  with  the  directions  of  our  future  work  on  this  subject. 

2.  MATHEMATICAL  AND  NUMERICAL  MODEL 

2.1  Objective  of  the  Model 

Detailed  modeling  of  gas-phase  reactive  flows  is  based  on  a  generally  accepted  set  of  time- 
dependent  coupled  partial  differential  equations  maintaining  conservation  of  total  mass,  momentum, 
total  energy  and  individual  species  density.  These  equations  describe  the  convective  motion  of  the 
fluid,  the  chemical  reactions  among  the  constituent  species,  and  the  diffusive  transport  processes 
such  as  thermal  conduction  and  molecular  diffusion. 


2.2  Gas  Phase  Equations 


A  strong  conservation  form  of  the  two-dimensional,  unsteady,  compressible  Navier-Stokes  equa¬ 
tions,  used  to  describe  gas  phase  reactive  flow  systems  can  be  written  as  follows  [27],  [28],  [29]. 


dp  dpu  dpv 

=  ^evap,  (1) 


-I-  P)  d(puv)  _  dr^x  ,  dTxy  ,  A 
dt  dx  dy  dx  dy 


,  djpuv)  d[pv^  -t-  P)  _  dTyx  dTyy  . 

dt  dx  ^  dy  dx  dy  + 


(3) 


^jp^)  ,  d{{pE+P)u)  d{{pE+  P)v)  _  djTxxu)  djrxyu)  d(ryxv)  d(Tyyv) 

dt  dx  dy  dx  dy  dx  dy 


dOx  dOy 

d  P9‘^  "t"  Qchem  "b  Qrad  "b  Q 


evapi 


dp^  ^  d{pku)  I  d{pkv)  _  d{pkUk)  d{pkVk)  ,  . 

dt  dx  dy  dx  dy  +  + 


A:  =  l,2,...,Ar. 


(4) 

(5) 


In  these  equations  x  and  y  denote  the  independent  spatial  coordinate  and  t  denotes  the  temporal 
coordinate;  p  the  mass  density;  pk  the  density  of  the  kth  species;  u  and  v  are  the  x  and  y  components 
of  the  fluid  bulk  velocity;  P ,  the  pressure;  E,  the  total  energy  of  the  fluid  per  unit  mass;  g,  the 
gravitational  body  force  per  unit  mass;  Uk  and  Vk  are  the  x  and  y  components  of  the  diffusion 
velocity  for  the  kth  species;  Txxi  Txyi  Tyy  are  the  components  of  the  stress  tensor  for  newtonian  fluid 
in  rectangular  coordinate;  qx  and  qy  are  the  x  and  y  components  of  the  heat-flux  vector,  mgvap  is 
the  source  term  arising  in  the  mass  conservation  equation  due  to  evaporation  of  the  water  droplets, 
resulting  in  the  production  of  water  vapor.  The  corresponding  term  in  the  H2O  gas  phase  species 
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conservation  equation  is  represented  by  mk,evap-  is  the  rate  of  production  of  the  species  due 
to  chemical  reactions.  Similarly  Qchem  is  an  exothermic  source  term  in  the  energy  equation  arising 
due  to  the  oxidation  of  fuel  molecules,  whereas  Qeyap  represent  absorption  of  energy  due  to  drop 
evaporation;  Qrad  is  the  radiative  heat  loss  term.  Sx  and  Sy  are  the  source  terms  arising  in  the  x 
and  y  momentum  equations  representing  the  cumulative  drag  force  exerted  by  the  various  droplet 
sections  on  the  gas  phase. 


The  viscous  stress  terms  are  related  to  the  rate  of  strain  through  the  stress-strain  relationship, 


4  du  2  dv'\ 


'^XX  —  M 

Zdx  3%J 

'du  dv' 

Tjjy  - 

.dy  dx. 

'4  dv  2  du' 

,Zdy  Zdx. 

(6) 


where  //  is  the  dynamic  viscosity.  The  components  of  the  heat  flux  vector  contain  contribution  due 
to  the  heat  conduction  and  the  interdiflfusion  process,  and  can  be  expressed  as 

qx  =  +  '^{hkPkUk)i 

k=i 

q^  =  -X^  +  'p{hkPkVk),  (7) 


where  A  is  the  thermal  conductivity  of  the  fluid  mixture  and  hk  is  the  enthalpy  per  unit  mass  for 
the  kih  species. 


The  total  energy  per  unit  mass  E  is  defined  as  the  sum  of  internal  energy  e  and  the  kinetic 
energy  as  follows 

E  =  e  +  ^{u'^  +  (8) 

The  internal  energy  is  related  to  the  thermodynamic  pressure  P  through  the  ideal  gas  equation  of 
state 


N 


p  =  r,tY. 


fc=l 


pk 

Wk’ 


and  the  caloric  equation  of  state 


e 


N 


k^i  p 


hk  = 


(9) 


(10) 

(11) 


In  these  equations  Ro  is  universal  gas  constant;  Wk,  the  molecular  weight  of  the  species  and 
Ih®  heats  of  formation  for  the  various  species  at  the  reference  temperature  T°. 
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2.3  Thermodynamics 


The  diffusion  velocities  Uk  and  14  for  the  species  are  computed  using  Pick’s  Law,  without 
solving  the  complete  matrix  problem  associated  with  a  multi-component  system. 

^k  =  -^Vkm^  k=l,2,...,N, 

1  dX 

~  ^ Qy  k  —  1,  2,  .  .  iV.  (12) 

Here,  Vkm  is  the  diffusion  coefficient  for  species  k  diffusing  into  the  mixture  and  Xk  is  the  mole 
fraction  of  the  k*^  species.  The  resulting  diffusion  velocities  are  corrected  to  satisfy  the  requirement 
that  the  diffusion  velocities  do  not  introduce  any  net  bulk  velocity  to  the  fluid.  The  diffusion 
velocities  are  computed  subject  to  the  constraint  equations 

N 

^PkUk  =  0, 
k=zl 
N 

PkVk  =  0,  (13) 

A;=l 


The  mole  fractions  can  be  related  to  the  mass  fraction  14  through  the  auxiliary  relation 

YkIWk 


Xk  = 


EtiYk/Wk' 


A:  =  l,2,...,iV. 


(14) 


The  mixture  diffusion  coefficient  'Dkm  of  the  species  is  related  to  the  binary  diffusion 
coefficient  V^i  of  species  k  diffusing  into  species  /.  The  binary  diffusion  coefficients  are  calculated 
from  kinetic  theory  and  are  in  the  following  form 

Vki  =  AkiT^^^jn,  A;,/=  l,2,...,iV,  (15) 

where  Aki  and  Bki  depend  on  species  k  and  /.  Values  for  Aki  and  Bki  have  been  tabulated  by 
Kailasanath  et  al.  [30].  The  diffusion  coefficient  of  species  A:  in  a  mixture  of  N  species  is  calculated 
according  to 


^km 


l-Yk 


Xi 

\jl^k  I 


-,/?  =  !, 2,..., AT. 


(16) 


The  mixture  thermal  conductivity  A  is  obtained  by  combining  the  thermal  conductivities  of 
the  individual  gases  Xk  that  are  present  in  the  mixture.  The  species  thermal  conductivity,  A^,  are 
estimated  theoretically  from  kinetic  theory  over  the  temperature  range  300/4"  to  3300/4",  and  these 
values  were  fit  to  a  third-order  polynomial.  The  mixture  thermal  conductivity  is  then  calculated 
using  the  expression  from  Kee  et  al.  [31] 


A 


1 

2 


1 

S  ^kXk  +  — 

jk=:l  2^k=l  Afc 


(17) 
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The  dynamic  viscosity  of  a  single  species  pik  can  be  obtained  from  kinetic  theory  and  is  fitted  over 
a  suitable  temperature  range  using  a  third  order  polynomial.  The  mixture  dynamic  viscosity  is 
calculated  using  the  expression  by  Wilke  [32], 


N 

/•  =  E 


XkFk 


where, 


m 

Wi 


(18) 


(19) 


2.4  Kinetics 

Ideally,  we  would  like  to  simulate  the  chemical  reactions  by  including  a  detailed  set  of  elemen¬ 
tary  reactions  to  describe  the  production  of  the  individual  species  and  the  energy  release  in  the 
flame.  However,  the  cost  of  computer  time  and  memory  required  to  track  the  individual  species 
makes  this  prohibitive  for  problems  in  which  parametric  studies  are  planned.  Instead,  the  chemical 
reaction  and  energy-release  process  for  methane-air  combustion  is  described  phenomenologically 
based  on  a  single  step  reaction, 

CH4  +  2O2  CO2  +  2H2O,  (20) 

using  a  finite-rate,  quasi-global  Arrhenius  expression.  Westbrook  et  al.  [33]  has  developed  simpli¬ 
fied  reaction  mechanisms  for  the  oxidation  of  hydrocarbons  using  a  laminar  flame  model  and  has 
prescribed  an  overall  rate  of  consumption  of  methane  as 

=  -1.3E8expi-48m/RT)[CH4]-'^-^[02Y-^.  (21) 

at 

The  above  expression  is  used  to  compute  the  depletion  of  methane.  The  stoichiometric  coefficients 
provided  in  Equation  (20)  is  then  used  to  compute  the  corresponding  changes  in  the  concentration 
of  oxygen,  carbon  dioxide,  water  vapor  and  nitrogen.  The  heat  release  term  Qchem  appearing  in 
the  energy  equation  4  is  determined  from 

Q,hem  =  (22) 

where  AHc  is  the  heat  of  combustion  for  the  single  step  global  chemical  reaction  per  unit  mass 
of  fuel.  The  heat  of  combustion  is  computed  using  the  stoichiometric  coefficients  and  the  heat  of 
formation  of  the  various  species  as  follows, 

AH,  =  j2i^'k-ui)Ahlk,  (23) 

k=l 


where,  and  are  the  stoichiometric  coefficients  of  the  species  as  a  reactant  and  product 
respectively  and  AHj  f.  is  the  heat  of  formation  of  the  species. 
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2.5  Water  mist  model 

In  this  section  we  describe  the  approach  that  we  have  used  in  formulating  the  system  of  dif¬ 
ferential  equations  that  govern  the  spray  and  gas  behavior  in  a  coupled  manner.  We  will  consider 
a  two-continuum  formulation,  wherein  the  gas  properties  and  the  droplet  properties  are  each  de¬ 
scribed  by  equations  in  the  Eulerian  form.  This  Eulerian-Eulerian  form  is  useful  when  resolution  is 
desired  only  on  a  scale  larger  than  the  average  distance  between  the  droplets.  In  this  approach  the 
droplet  properties  are  treated  as  if  they  were  continuous  in  the  domain  as  the  gaseous  properties 
[34],  [35].  The  droplet  property  at  a  point  in  space  and  time  represents  the  average  value  over  many 
droplets  in  the  neighborhood  of  that  point. 

The  size  distribution  of  the  spray  droplets  can  be  described  by  the  concentration  of  discrete 
droplets  of  various  sizes  per  unit  volume  of  fluid  ni{x,y,t)  as  a  function  of  spatial  coordinates  x, 
y  and  of  time  t  where  i  =  1,2,....  The  magnitude  of  the  integer  i  represents  the  total  number 
of  discrete  droplet  sizes.  The  vaporization  process  [36],  [37]  can  be  described  by  a  set  of  coupled 
differential  equations  for  the  concentration  of  the  discrete  droplet  sizes 

dui  dui  dui  ^  ^ 

dt  ^  '^^'^~dy  ~  >  *=1,2,...,  (24) 

where  and  u/,,-  are  the  x  and  y  components  of  the  velocity  vector  of  the  droplet  size  and  Ei 
is  the  frequency  of  the  molecule  evaporation.  This  frequency  highly  depends  on  the  surface  area  of 
the  droplet  and  on  the  temperature  of  the  ambient  gas  in  addition  to  other  flow  parameters.  Since 
the  total  number  of  droplet  sizes  needed  to  simulate  actual  fuel  sprays  can  be  immense,  sectional 
conservation  equations  have  been  developed.  This  method  is  based  on  dividing  the  droplet  size 
domain  into  sections  and  dealing  only  with  one  integral  quantity  in  each  section.  This  sectional 
representation  has  the  advantage  that  the  integral  quantity  is  conserved  within  the  computational 
domain  and  the  number  of  conservation  equations  is  substantially  reduced  so  as  to  equal  the  number 
of  sections. 


We  divide  the  entire  droplet  size  domain  into  M  arbitrary  sections,  and  define  Qj  to  be  an 
integral  quantity  of  the  spray  with  the  section.  Thus 


i  =  l,2,...,M, 


(25) 


where  n(u,  t)  is  the  number  concentration  function  and  Vj^i  and  Vj  denote  the  volumes  of  the 
smallest  and  largest  droplets,  respectively,  in  section  j.  Using  equation  25  in  equation  24  the 
sectional  conservation  equation  can  be  rewritten  as 


dQj  dQi  dQi  _ 

"df  ^  ,  i  =  1, 2, . . . ,  M, 

where  the  sectional  coeflficients  Bjj+i  and  Cj  are  described  in  the  following  forms: 


(26) 


B 


’jj+i 


_  f  Vkj  \  E{vj,T) 


Vk^+l)  Vj+i 


-  Vi 


(27) 


Cj  =  (  ^  + - - -  r  -E{v,T)dv, 

\Vkj-i+lJ  Vj-Vj-i  Vj  -  Vj^i  V 


Vkj-1+1 


(28) 


Water  Mist  Fire  Suppression 


9 


where  E{v,T)  is  the  volume  rate  of  vaporization  of  a  droplet  of  volume  v. 

The  mass  or  volume  rate  of  vaporization  of  each  single  droplet  is  proportional  to  the  droplet 
diameter  d  (the  d?  law).  Thus 


Eiv,T)  =  {7r/4)EiT)d,  (29) 

where  E{T)  is  the  vaporization  or  burning  rate  coefficient,  which  depends  on  the  temperature  of  the 
surrounding  gas,  the  diffusivity  and  other  properties  of  the  droplet  and  its  surroundings.  For  the 
present  calculations,  the  burning  rate  coefficient  value  of  7.6EZiJ,m^/s  was  assumed  for  fine  water 
mist  particles  [38],  [27].  Assuming  a  d^  law  to  re-express  the  sectional  vaporization  coefficients  in 
terms  of  droplet  diameters  and  assuming  a  continuous  division  into  sections  the  coefficients  can  be 
rewritten  as 


=  -E{T) 


Q  =  ^EiT) 


dL,j+i 


idu,j+i  -  d%j+ij 


d-Lj 


y  =  1,2,...,M-  1, 
3(d[7j  -  dij) 


ldu,j  -  dlj 


du,j-dl,j  J 


j  =  1,2,  ...,M. 


(30) 


Each  of  the  droplet  sections  is  assumed  to  have  its  own  unique  velocity  different  from  that  of 
the  gas  phase.  Momentum  conservation  equations  are  formulated  for  each  droplet  section  and  are 
coupled  to  those  of  the  gas  phase  through  the  phase  interaction  terms  (drag  terms).  The  droplet 
drag  equations  are  given  by 


duij 

dt 


+  utj 


duij 

dx 


dvi,i 

dt 


+ 


dx 


+  ‘^1,3 


+  VI, 3 


dui,j 

dy 

dvi,3 

dy 


—  Ex,j )  j  —  1)2,...,  Ml 

~  Ey,j  ~  y,  J  =  1)  2,  •  ■  • ,  M, 


(31) 


where  Fxj  and  Fyj  are  the  x  and  y  components  of  the  drag  force  acting  on  the  sectional  velocity 
modeled  based  on  Stokes  Law.  The  gas  phase  momentum  source  terms  Sx  and  Sy  are  equal  and 
opposite  to  the  sectional  density  weighted  sum  of  Fxj  and  Ey,3  respectively. 


2.6  Numerical  Algorithm 


The  governing  equations  are  rewritten  in  terms  of  finite-volume  approximations  on  an  Eulerian 
mesh  and  solved  numerically  for  specific  boundary  and  initial  conditions.  A  complete  solution  to 
these  governing  equations  require  solving  the  terms  for  each  of  the  individual  processes,  as  well 
as  accounting  for  the  interaction  among  the  processes.  The  solution  approach  consists  of  separate 
algorithms  for  each  of  the  individual  processes,  which  are  then  coupled  together  by  the  method 
of  time-step  splitting.  The  algorithms  for  convection,  thermal  conduction,  molecular  diffusion, 
viscosity  and  the  coupling  of  the  individual  processes  have  been  previously  discussed  in  detail  [39], 
[40]. 


The  fluid  convection  is  solved  with  a  high-order  implicit  algorithm.  Barely  Implicit  Correction  to 
the  Flux-Corrected  Transport  (BIC-FCT)  [39],  that  was  developed  to  solve  the  convection  equations 
for  low-velocity  flows.  The  Flux-Corrected  Transport  (FCT)  [41]  algorithm  itself  is  an  explicit. 
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finite-difference  algorithm  that  is  constructed  to  have  fourth-order  phase  accuracy.  Through  a  two- 
step  predictor-corrector  algorithm,  FCT  ensures  that  all  conserved  quantities  remain  monotone  and 
positive.  The  FCT  procedure  is  to  first  modify  the  properties  of  a  high-order  algorithm  by  adding 
diffusion  during  a  convection  step  and  then  to  subtract  out  the  diffusion  in  an  antidiffusion  phase. 
In  addition,  fluxes  are  limited  to  ensure  that  no  new  unphysical  maxima  or  minima  are  added 
during  the  convection  process.  However,  because  FCT  is  an  explicit  algorithm,  the  numerical  time- 
step  required  for  accuracy  and  stability  is  limited  by  the  velocity  of  sound  according  to  the  CFL 
condition.  To  filter  out  the  sound  waves  from  the  convection  equations  and  therefore  remove  the 
sound-speed  (Courant)  limitation  on  the  time-step,  the  convection  equations  are  usually  solved 
implicitly.  The  BIC-FCT  algorithm  was  developed  so  that  the  time  step  is  limited  by  the  fluid 
velocity  and  not  the  sound  speed.  This  implementation  is  particularly  useful  for  slowly  evolving 
flows  because  one  BIC-FCT  time-step  requires  the  same  amount  of  computer  time  as  one  regular 
FCT  explicit  time  step,  but  the  size  of  the  timestep  might  be  a  factor  of  50-100  times  greater. 
The  reaction  rate  equations  involving  the  arrhenius  terms  are  solved  using  an  implicit  trapezoidal 
scheme  with  time  step  splitting  over  the  largest  fluid  dynamic  time  step. 


A  specific  solution  of  the  reactive  flow  equations  is  determined  by  the  initial  conditions  and  the 
boundary  conditions  that  describe  the  geometry  of  the  system  and  exchange  of  mass,  momentum 
and  energy  occurring  between  the  system  and  the  rest  of  the  physical  world.  The  total  density, 
u  and  V  momentum  and  the  various  species  densities  are  prescribed  at  the  inflow  boundary.  At 
the  outflow  boundary,  the  normal  gradients  of  total  density,  species  densities  and  momentum  was 
assumed  to  be  zero.  The  pressure  at  large  distances  from  the  burner  surface  is  assumed  to  be  equal 
to  the  ambient  pressure.  A  symmetric  boundary  condition  was  employed  at  the  center  line  of  the 
computational  domain  by  use  of  anti-symmetric  reflection  of  tangential  velocity  u,  and  symmetric 
reflection  of  all  other  variables,  A  slip  wall  boundary  condition  is  employed  at  the  lateral  boundary 
by  assuming  that  the  flux  of  all  transported  properties  across  the  wall  is  zero.  For  a  diflFusion  flame 
computation,  the  initial  conditions  at  each  point  in  the  domain  were  the  same  as  inlet  conditions. 
The  spray  injection  densities  and  the  two  components  of  the  velocity  vector  are  specified  at  the 
inflow  boundary.  Again  symmetric  boundary  conditions  are  employed  for  the  liquid  phase  at  the 
center  line  and  slip  wall  conditions  are  employed  at  the  edge  of  the  burner. 


The  computational  domain  consists  of  a  stretched  96  X  128  grid,  concentrated  near  the  contact 
point  of  the  fuel  and  the  oxidizer.  The  smallest  grid  cell  is  a  0,5mm  x  0.5mm  square  cell.  Numerical 
simulations  were  performed  by  halving  the  grid  size  to  check  for  accuracy  and  consistency  of  the 
algorithms.  The  algorithms  have  been  optimized  for  efficient  performance  on  vector  machines. 
Typical  computational  time  on  a  C-90  single  processor  machine  was  approximately  4  hours  to 
reach  a  steady  state  solution. 


3*  RESULTS  AND  DISCUSSION 


The  governing  equations  described  earlier  were  solved  using  the  numerical  models  briefly  dis¬ 
cussed  above,  to  study  methane-air  diffusion  flames  stabilized  above  a  Wolfhard  Parker  burner. 
Fine  water  mist  is  introduced  along  with  the  co-flowing  air  to  study  the  interaction  and  subse¬ 
quent  suppression  of  these  flame  with  water  mist.  The  numerical  results  have  been  compared  with 
experimental  data  obtained  using  a  similar  geometry. 
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3.1  Comparison  to  Experimental  Data 

The  key  component  of  the  experimental  setup  is  a  modified  Wolfhard-Parker  burner  which 
is  shown  schematically  in  Figure  1.  The  burner  slot  is  75  mm  long,  10  mm  wide  and  150  mm 
deep.  It  has  two  identical  oxidizer  channels  82  mm  long,  35  mm  wide  and  150  mm  deep  on  each 
side.  Beside  each  oxidizer  channel  is  the  mist  generation  chamber  where  mist  is  produced  with 
commercial  low  flow  Delavan  nozzles.  A  fraction  of  the  mist  generated  in  this  chamber  is  entrained 
into  the  air  stream.  The  slot  opening  is  adjustable  to  control  the  quantity  of  mist  entrained  into 
the  air  stream.  Mist  flow  rate  is  determined  by  covering  the  channel  for  a  given  time  with  a  very 
fine  screen.  The  mass  of  mist  collected  during  this  period  gives  a  measure  of  the  mist  flow  rate  out 
of  the  air  channel.  To  introduce  steam  along  with  the  air  co-flow,  an  electrically  heated  fine  screen 
is  used  to  collect  and  evaporate  the  mist  droplets  before  they  exit  from  the  air  channel.  The  fuel 
slot  is  filled  with  fine  clean  sand  to  produce  a  near  plug  flow  condition.  Similarly,  three  layers  of 
fine  screen  are  used  to  produce  a  uniform  air  velocity  profile  in  each  air  channel.  In  the  current 
experiments,  the  average  exit  velocity  of  the  fuel  is  2.81cm/sec  and  the  average  exit  velocity  of  the 
air  is  18.1cm/sec.  A  quartz  glass  shield,  47  cm  long  and  27  cm  high  sits  on  the  burner  on  either  side 
to  protect  the  flame  from  any  external  disturbances.  The  flame  temperature  is  measured  with  a  50 
fim  diameter  fine  platinum/platinum  -  13%  -  rhodium  thermocouple  which  is  computer  controlled 
by  a  Newport  3D  positioning  instrument  The  droplet  diameter  is  characterized  using  the  Malvern 
particle  analyser.  Infrared  images  of  the  flame  are  also  obtained  using  the  Agema  Thermovision  870 
scanner,  calibrated  with  a  blackbody  source  and  with  the  thermocouple  temperature  measurements. 

Figure  2  shows  comparison  of  numerically  computed  and  experimentally  measured  temperature 
profiles  for  the  the  base  case  methane  -  air  diffusion  flames.  The  comparison  has  been  made  at  a 
height  of  5,  15,  30  and  55  mm  above  the  burner  surface.  At  each  height  the  numerically  computed 
temperature  profile  compares  favorably  with  the  thermocouple  temperature  measurements  [43]. 
The  numerically  computed  values  are  very  sensitive  to  the  boundary  conditions  that  are  imposed 
at  the  burner  surface.  The  small  differences  in  temperature  profiles  between  the  numerical  and 
experimental  results  are  due  to  the  idealization  of  the  boundary  conditions  observed  in  the  experi¬ 
ments.  It  is  not  possible  to  mathematically  characterize  the  precise  boundary  conditions  present  in 
the  experiments.  Therefore,  experimentally  determined  values  were  used  as  boundary  conditions 
in  the  numerical  solution  procedure.  Once  set,  the  same  boundary  conditions  have  been  used  for 
all  the  cases  discussed  in  this  paper. 

Next,  the  effect  of  nitrogen  dilution  on  methane  -  air  diffusion  flame  was  studied.  Nitrogen 
was  introduced  along  with  the  air  co-flow  resulting  in  different  levels  of  nitrogen.  Figure  3  shows 
comparison  of  numerical  and  experimental  results  for  CH^  and  18.2%02, 81.8%iV2  mole  fraction 
diffusion  flame.  Numerically  computed  temperature  profiles  are  compared  with  experimental  re¬ 
sults  at  a  height  of  6,  40,  46  and  50  mm  above  the  burner  surface.  Again  the  agreement  between 
the  computed  and  experimentally  measured  values  is  very  good.  Additional  cases  with  different 
amounts  of  nitrogen  dilution  (Figure  4  for  the  case  of  16.0%O2,  83.5%iV2  mole  fraction  )  have  also 
been  simulated  and  have  been  compared  with  experimental  results.  As  nitrogen  dilution  increases 
the  diffusion  flame  appears  to  become  broader  and  taller  due  to  a  shift  in  the  location  of  the  stoi¬ 
chiometric  surface  towards  the  oxidizer  jet.  As  expected,  the  maximum  flame  temperature  reduces 
with  increasing  levels  of  nitrogen  dilution  until  a  critical  dilution  level  at  which  the  flame  extin¬ 
guishes.  This  effect  is  similar  to  the  “oxygen  displacement  effect”  observed  during  the  interaction 
of  water  mist  with  diffusion  flame  (discussed  later  in  the  text).  The  exact  level  of  nitrogen  dilution 
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for  extinction  can  be  predicted  more  accurately  by  using  detailed  finite  rate  kinetics  for  modeling 
the  oxidation  of  the  fuel  molecule.  Due  to  the  very  high  computational  costs  cissociated  with  using 
finite  rate  kinetics,  it  is  not  possible  to  perform  the  large  number  of  parametric  studies  required 
for  optimizing  the  water  mist  spray  characteristics.  The  use  of  single  step  chemistry  is  shown  to 
reproduce  the  correct  trends  as  observed  in  the  experiments,  while  keeping  the  computational  costs 
low. 

3.2  Base  Case  Diffusion  Flame  Results 

One  of  the  axivantages  of  numerical  simulations  is  that  a  detailed  spatial  characterization  of  the 
flame  and  flow  field  can  be  obtained  with  little  additional  cost  or  effort.  In  addition,  quantities  such 
as  heat  release  rate  can  be  predicted  which  are  difficult  to  obtain  through  experiments.  Such  detailed 
numerical  diagnostics  will  be  used  to  understand  the  interaction  between  water  mist  and  the  flame. 
Figure  5  shows  contour  plots  of  temperature  {K)  and  heat  release  rate  {J oules / mtr^ / sec)  above  a 
two-dimensional  methane-air  diffusion  flame  burner.  These  numerical  simulations  were  performed 
for  the  burner  geometry  and  flow  conditions  described  earlier.  The  temperature  contours  indicate 
the  presence  of  a  very  thin  flame  sheet  which  closes  above  the  fuel  duct.  The  height  of  the  flame 
is  roughly  7cm.  The  maximum  flame  temperature  is  1970A',  which  agrees  well  with  experimental 
results  [43],  [44].  The  heat  release  rate  profile  indicates  the  region  of  rapid  exothermic  chemical 
activity.  This  profile  identifies  the  location  of  the  flame  sheet  where  the  fuel  and  the  oxidizer 
species  come  together  in  stoichiometric  proportions,  react  exothermically  and  produce  product 
species.  Figure  6  shows  contour  plots  of  methane  {CH4)  and  oxygen  (O2)  density  {kg/mtr^)  above 
the  burner  surface.  Due  to  the  fuel  lean  conditions  in  which  these  simulations  were  performed, 
methane  is  completely  consumed.  The  unused  oxygen  flows  out  of  the  computational  domain 
along  with  other  product  species  and  diluent  gases.  Figure  7  shows  density  profiles  of  the  product 
species,  carbon-dioxide  (CO2)  and  water  vapor  [H2O).  Figure  8  shows  the  velocity  vectors  and 
streamline  pattern.  The  velocity  vectors  have  been  color  coded  with  red  showing  the  highest  gas 
speed  and  blue  showing  the  lowest  gas  velocity.  The  direction  of  the  arrows  shows  the  direction  of 
the  velocity  vector.  The  streamline  pattern  clearly  indicates  the  entrainment  of  oxidizer  along  with 
the  suppressants  into  the  diffusion  flame.  Both  the  streamlines  and  the  velocity  vectors  bend  into 
the  flame  due  to  the  entrainment  of  air  into  the  diffusion  flame.  When  similar  calculations  were 
performed  with  the  buoyancy  terms  set  to  zero,  no  such  entrainment  or  bending  of  the  streamlines 
into  the  flame  Wcis  observed.  Infact,  the  streamlines  under  zero  gravity  conditions  tend  to  diverge 
away  from  the  flame  due  to  expansion  of  the  hot  gases. 

3.3  Interaction  of  Water  Mist  with  Flame 

The  primary  goal  of  the  present  research  is  to  study  the  interaction  of  water  mist  with  diffusion 
flames.  With  the  confidence  gained  by  the  comparison  with  experimental  data,  a  series  of  numerical 
simulations  were  performed  to  determine  the  underlying  processes  involved  in  suppression  of  flames 
with  water  mist.  Water  mist  was  introduced  along  with  the  air  co-flow.  The  spray  velocity  at  the 
burner  exit  was  set  at  25cmf  sec.  Figure  9  shows  temperature  and  heat  release  rate  profiles  for  50/x 
droplets  injected  uniformly  over  the  length  of  the  air  channel.  The  spray  injection  density  was  held 
fixed  at  2000drops/cTn^.  The  maximum  flame  temperature  reduces  from  1970A'  in  the  base  case  to 
1795A'  when  subjected  to  50ju  water  droplets.  The  predicted  flame  height  increases  from  7cm  for  the 
base  case  to  9cm  when  subjected  to  fine  water  mist.  The  temperature  and  heat  release  rate  profiles 
indicate  an  increase  in  flame  height  and  flame  suppression  due  to  the  sub-extinguishing  application 
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rates  of  water  mist.  The  heat  release  rate  profile  also  indicates  the  presence  of  an  endothermic 
heat  release  rate  profile  (negative  value)  showing  the  approximate  location  of  evaporation  of  the 
droplets.  When  the  droplets  are  introduced  parallel  to  the  air  co-flow,  the  evaporation  time  at 
the  flame  sheet  was  much  smaller  than  the  rate  of  entrainment.  For  the  conditions  that  have 
been  tested,  we  did  not  observe  any  droplets  that  could  survive  the  very  high  temperatures  at  the 
flame  and  penetrate  through  the  diffusion  flame.  Figure  10  shows  similar  contours  of  temperature 
and  heat  release  under  the  application  of  150/Li  water  droplets  and  a  spray  injection  density  of 
lOOdrops / cm^.  Flames  subjected  to  water  mist  are  taller  and  broader  than  the  base  case  methane- 
air  diffusion  flames.  The  flame  height  increases  by  approximately  2cm  for  this  case.  The  maximum 
flame  temperature  reduces  from  1970*'  in  the  base  case  to  1905*"  when  subjected  to  150/i  water 
droplets.  The  endothermic  values  (-ve  numbers)  in  the  heat  release  rate  contours  are  approximately 
two  orders  of  magnitude  smaller  than  the  exothermic  values.  However  the  area  occupied  by  the 
endothermic  contours  is  much  larger  than  that  of  the  exothermic  regions.  Although  the  evaporation 
time  is  proportional  to  the  square  of  the  droplet  diameter,  the  characteristic  time  for  decrease  of  the 
relative  velocity  due  to  the  drag  forces  is  also  proportional  to  the  square  of  the  droplet  diameter. 
As  a  result  with  increasing  droplet  diameter  we  do  not  expect  the  droplets  to  penetrate  the  high 
temperatures  of  the  flame  sheet  as  long  as  they  are  introduced  in  a  co-flow  manner. 


Numerical  simulations  have  been  performed  to  ciscertain  the  effect  of  changing  the  injection 
spray  density  on  methane- air  diffusion  flames.  Figure  11  shows  temperature  contours  for  various 
spray  injection  densities  of  1000,  2000,  3000  and  4000  drops/cm^  of  50/x  water  mist  droplets.  As 
in  the  previous  simulations  the  spray  is  injected  along  with  the  air  co-flow  perpendicular  to  the 
burner  surface.  As  the  spray  injection  density  increases  more  and  more  water  droplets  entrain 
into  the  diffusion  flame  and  this  results  in  larger  suppression  as  shown  by  the  reduction  in  flame 
temperature.  We  also  observe  that  as  the  spray  injection  density  increases,  the  flame  becomes 
taller  and  broader  (fatter),  resulting  in  larger  flame  surface  area.  This  result  is  similar  to  the 
one  observed  under  nitrogen  and  steam  dilution  of  the  diffusion  flame.  The  flame  height  increases 
from  7cm  for  the  base  cases  flame  to  about  9cm  for  the  2000drops/cm®  case  and  to  11cm  for  the 
injection  of  4000drops/cm^  into  the  co-flow.  The  initial  injection  velocity  of  the  droplets  was  held 
fixed  at  2bcm.fsec  for  all  spray  densities.  The  droplets  are  divided  into  five  sections  (as  discussed 
in  section  dealing  with  mathematical  formulation)  as  follows  :  0  —  lOp,  10  —  20p,  20  —  30/i,  30  —  40^ 
and  40  -  50//.  Droplet  section  densities  for  three  of  the  five  sections  have  been  shown  in  Figure 
12.  These  contours  show  the  location  of  maximum  density  and  position  where  these  droplets 
evaporate  to  form  smaller  and  smaller  droplets  (sections).  The  injected  spray  is  a  monodisperse 
spray  and  consists  entirely  of  50//  droplets.  These  droplets  convect  along  with  the  air  co-flow. 
Some  of  these  droplets  evaporate  and  form  smaller  droplet  sections  while  others  appear  to  convect 
out  of  the  computational  domain.  The  droplets  that  entrain  into  the  diffusion  flame  evaporate 
and  form  smaller  sections,  resulting  in  a  very  large  density  of  the  smaller  sections  very  close  to 
the  diffusion  flame.  As  a  result,  the  profiles  show  a  very  distinct  cascading  effect  in  which  larger 
diameter  droplets  show  maximum  concentration  further  away  from  the  diffusion  flames,  whereas 
smaller  diameter  droplets  show  maximum  concentration  closer  to  the  diffusion  flame. 

Figure  13  shows  similar  droplet  section  density  contours  for  the  injection  of  150//  droplets. 
In  this  case  the  droplets  are  again  divided  into  five  sections  of  0  —  30//,  30  —  60//,  60  —  90//, 
90  -  120//  and  120  -  150//.  Again  the  injected  spray  is  a  monodisperse  spray  consisting  entirely  of 
150//  droplets.  Unlike  the  earlier  case  of  50//  droplets,  in  this  case  the  120  -  150//  section  advects 
with  the  gas  phase  only  for  a  very  small  distance  above  the  burner  surface.  Beyond  this  point 
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the  section  density  and  section  velocity  are  zero.  This  is  because  the  drag  force  acting  on  these 
droplets  is  not  able  to  counter  balance  the  weight  of  these  droplets  and  as  a  result  the  droplet 
velocity  continuously  reduces.  When  the  droplets  are  injected  at  higher  initial  inlet  velocity,  it  was 
observed  that  the  droplets  were  able  to  reach  a  higher  distance  above  the  burner  surface  due  to  the 
larger  initial  momentum.  The  smaller  sections  again  show  profiles  similar  to  those  observed  with 
the  50^  droplet  case. 

3.4  Droplet  Dynamics 

The  numerical  model  was  used  to  obtain  a  detailed  understanding  of  the  droplet  dynamics 
in  a  2-dimensional  flow  field  and  interaction  /  entrainment  of  mist  into  the  diffusion  flame.  As 
described  earlier  the  droplet  range  is  divided  into  5  sections.  The  movement  of  the  various  droplet 
sections  was  analyzed  by  the  velocity  profiles.  Velocity  profiles  (velocity  component  normal  to  the 
burner  surface)  of  the  various  sections  along  with  the  gcis  phase  profile  have  also  been  shown  as  a 
function  of  distance  above  the  burner  surface  in  figure  14.  By  comparing  these  profiles,  the  effect 
of  droplet  drag  and  buoyancy  forces  on  the  various  droplet  sections  can  be  deduced.  These  results 
show  that  smaller  diameter  droplets  (50fi)  quickly  follow  the  gas  phase  (small  characteristic  time 
for  decrease  of  relative  velocity)  and  are  entrained  into  the  flame  at  approximately  the  same  rate  as 
the  oxidizer.  However  larger  diameter  droplets  (150//)  tend  to  travel  with  their  injection  velocity 
and  exhibit  larger  characteristic  times  for  the  decrease  in  relative  velocity.  For  the  smaller  droplets 
the  drag  force  exerted  by  the  gas  on  these  droplets  is  able  to  counterbalance  the  droplet  weight 
and  these  droplets  therefore  entrain  into  the  flame  along  with  the  oxidizer.  The  larger  droplets  are 
not  able  to  counterbalance  their  weight  with  the  drag  force.  Due  to  the  varied  characteristic  times 
associated  with  the  different  diameter  droplets  it  is  found  that  the  effect  of  injection  characteristics 
on  flame  suppression  are  more  critical  for  the  larger  diameter  droplets.  The  smaller  diameter 
droplet  entrain  most  rapidly  into  the  flame  when  injected  in  a  co-flow  configuration,  whereas  larger 
diameter  droplets  entrain  at  a  slower  rate. 

3.5  Comparison  of  Suppressant  Results  to  Experimental  Data 

Figure  15  shows  the  comparison  between  numerically  computed  and  experimentally  measured 
temperature  profiles  in  a  methane-air  diffusion  flame  when  steam  (3%  mass  fraction)  is  added 
to  the  co-flow  air.  The  comparison  has  been  made  at  the  heights  of  10mm,  30mm  and  45mm 
above  the  burner.  Figure  16  shows  a  similar  comparison  when  water  mist  (3%  mass  fraction)  was 
added.  In  general  the  numerically  predicted  profiles  compare  well  with  the  experimentally  measured 
profiles.  However  the  comparisons  become  poorer  as  the  height  above  the  burner  increases.  An 
interesting  observation  from  figures  15  and  16  is  that  water  mist  suppresses  flame  temperature 
more  than  steam  does.  This  is  due  to  the  role  played  by  the  latent  heat  for  water  evaporation. 
The  mechanisms  of  flame  temperature  suppression  by  water  mist  are  discussed  in  more  detail  later 
in  this  section.  Experimentally  the  injection  of  water  mist  makes  the  flame  flicker.  This  results 
in  the  mixing  of  hot  and  cold  regions  near  the  flame  sheet  at  any  given  distance  above  the  burner 
surface.  Because  of  flickering,  it  is  difficult  to  define  a  steady  temperature  at  a  point  whenever 
mist  is  introduced.  Because  the  flame  moves,  the  thermocouple  encounters  various  regions  of  the 
flame  as  the  thermocouple  waits  to  make  a  measurement.  This  results  in  a  large  variation  of 
temperature  with  time.  The  effect  of  flickering  varies  with  location  in  the  flame  and  also  with  the 
mist  concentration.  At  lower  heights  and  near  the  center  of  the  flame,  it  is  not  as  severe,  but  it 
becomes  very  severe  near  the  edge  of  the  flame  as  the  plume  region  is  approached  (see  figure  17). 
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Such  changes  in  temperature  with  time  at  similar  locations  in  the  base  case  tests  and  in  tests  with 
nitrogen  are  insignificant.  At  each  point  where  temperature  is  measured,  maximum,  minimum  and 
mean  values  of  the  set  of  data  are  obtained.  The  experimental  data  in  figure  16  are  mean  time 
averaged  thermocouple  measurements.  Despite  the  uncertainty  resulting  from  flame  fluctuations, 
the  comparison  for  the  3%  mist  case  is  good. 

Figure  18  shows  comparison  of  numerically  predicted  and  experimentally  measured  temperature 
profiles  as  a  function  of  height  above  the  burner  surface.  The  numerical  result  is  the  maximum 
temperature  at  a  given  height  above  the  burner  surface,  whereas  the  experimental  results  include 
the  maximum,  mean  and  minimum  of  the  time  dependent  temperature  data.  At  each  height  above 
the  burner  surface,  the  time  dependent  profiles  are  experimentally  measured  at  various  points  across 
the  flame.  The  maximum,  mean  and  minimum  temperatures  shown  in  flgure  18  are  the  maximum 
values  at  a  given  height.  The  maximum  and  minimum  temperatures  at  a  given  height  thus  provide 
an  upper  and  lower  bound.  As  expected  the  numerical  results  fall  within  the  limits  predicted  by  the 
experimental  data.  The  maximum  possible  flame  suppression  was  also  computed  using  the  NASA 
thermochemical  equilibrium  code  [42].  As  expected,  the  numerically  observed  suppression  was  less 
than  the  maximum  possible  suppression  obtained  from  the  equilibrium  calculations. 

3.6  Mechanisms  of  Flame  Suppression 

There  are  four  basic  mechanisms  [43]  involved  in  the  extinction /suppression  of  flames  by  water 
mist.  These  include  thermal  cooling,  oxygen  displacement,  fuel  surface  cooling  and  radiation  atten¬ 
uation.  Since  the  focus  of  the  present  work  is  on  gas  phase  diffusion  flames,  the  effects  of  fuel  surface 
cooling  and  radiation  attenuation  have  been  neglected.  Numerical  simulations  were  performed  to 
examine  the  relative  contributions  of  thermal  cooling  and  oxygen  displacement  during  water  mist 
suppression  of  diflfusion  flames.  When  diffusion  flames  are  subjected  to  water-mist,  water  droplets 
evaporate  by  absorbing  energy  from  the  flame.  This  effect  of  the  latent  heat  of  evaporation  is 
included  in  the  “thermal  cooling”  effect.  Once  the  droplets  evaporate,  they  produce  water  vapor 
which  dilutes  the  oxidizer  or  the  surrounding  air  flow.  This  oxygen  displacement  further  suppresses 
the  flame  because  of  changes  in  rates  of  exothermic  chemical  reactions  because  of  changes  in  densi¬ 
ties  of  the  reactants.  The  water  vapor,  due  to  its  larger  heat  capacity  as  compared  to  air,  can  also 
change  the  specific  heat  of  the  the  reactant  gases  resulting  in  further  suppression  of  the  flame.  A 
more  drastic  effect  of  water  vapor  can  be  in  the  form  of  reduction  /  displacement  of  oxygen  to  a 
point  where  the  flame  extinguishes  due  to  oxygen  starvation.  This  scenario  is  more  relevant  to  fires 
that  occupy  a  large  part  of  their  enclosures.  All  these  effects,  other  than  that  associated  directly 
with  the  latent  heat  of  evaporation  are  included  in  the  “oxygen  displacement  effect” . 

As  described  earlier,  numerical  simulations  were  performed  with  50/r  droplets  injected 
along  with  the  air  co-flow  at  a  velocity  of  2bcmfsec  and  an  initial  spray  injection  density  of 
2000drops/c'm^.  These  results  were  shown  in  Figure  9.  This  numerical  simulation  was  repeated 
with  zero  value  for  the  heat  of  evaporation  to  examine  the  relative  contribution  of  the  thermal  cool¬ 
ing  and  oxygen  displacement  to  flame  suppression.  Figure  19  shows  temperature  profiles  for  a)  the 
base  case  methane-air  diffusion  flame,  b)  flame  subjected  to  water  mist  droplets  and  c)  flame  sub¬ 
jected  to  water  mist  droplets  but  with  latent  heat  of  evaporation  set  to  zero.  As  the  water  droplets 
are  introduced  into  the  flame  (Case  b),  the  flame  appears  to  cool  down  significantly  as  compared 
to  the  base  case  flame  (Case  a).  The  flame  also  becomes  taller  and  broader,  when  subjected  to  fine 
water  mist.  When  the  latent  heat  of  evaporation  is  set  to  zero  (Case  c),  the  flame  height  was  found 
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to  be  comparable  with  that  of  the  water  mist  case  (Case  b).  Figure  20  shows  the  temperature 
profiles  for  these  three  cases  at  four  heights  (25,  50,  75  and  IOOtthti)  above  the  burner  surface. 
At  the  heights  of  25mm  and  50mm  above  the  burner  surface,  the  base  case  temperature  profile 
almost  overlaps  the  temperature  profile  for  the  case  of  water  mist  with  zero  heat  of  evaporation. 
The  profile  for  water  mist  with  normal  values  of  heat  of  evaporation  shows  significant  suppression 
in  flame  temperature.  This  implies  that  at  heights  of  25mm  and  40mm  above  the  surface  the 
principal  effect  of  water  mist  is  through  the  latent  heat  of  evaporation  and  that  oxygen  dilution  has 
a  negligible  effect  on  flame  suppression.  As  one  moves  further  above  the  burner  surface,  the  effect 
of  oxygen  displacement  increases.  At  the  heights  of  75mm  and  100mm  above  the  burner  surface, 
the  numerical  results  indicate  that  oxygen  displacement  is  equally  important  for  flame  suppression 
as  is  thermal  cooling  [43],  [44].  Our  results  thus  indicate  that  the  effect  of  oxygen  dilution  is  more 
prominent  in  the  upper  part  of  the  flame,  whereas,  endothermic  heat  exchange  with  the  droplets 
play  a  critical  role  in  the  lower  half  of  the  diffusion  flame.  Figure  21  shows  temperature  profiles 
at  four  different  heights  above  the  burner  surface  for  the  three  cases  discussed  in  figure  19.  The 
droplet  diameter  in  this  case  was  150^  as  opposed  to  50p  for  Figure  20.  Again  the  results  indicate 
that  the  relative  contribution  of  the  various  suppression  mechanism  varies  with  height  above  the 
burner  surface  as  well  as  with  droplet  diameter  and  injection  velocity  of  the  water  droplets. 

3.7  Effect  of  Droplet  Diameter,  Injection  Density  and  Velocity 

A  series  of  parametric  studies  were  performed  to  estimate  the  effect  of  droplet  number  density 
and  droplet  mass  density  on  flame  suppression  and  flame  extinction.  Figure  22  shows  changes  in 
peak  temperature  as  a  function  of  droplet  mass  density.  Numerical  simulations  were  performed 
with  50,  100  and  150 fj,  droplets  with  varying  levels  of  injection  spray  densities.  The  injection 
velocity  of  the  droplets  was  unchanged  in  these  calculations.  The  peak  temperature  observed  for 
each  simulation  was  mapped  as  a  function  of  droplet  number  density  as  well  as  mass  density.  The 
profiles  of  peak  temperature  vs.  droplet  injection  mass  density  indicate  that  as  the  droplet  mass 
density  increases,  the  flame  temperature  reduces  until  a  certain  critical  mass  density  is  reached 
beyond  which  point  a  stable  flame  cannot  be  sustained.  This  trend  was  observed  for  all  water 
mist  droplet  diameters  studied  in  this  analysis.  It  is  also  observed  that  as  the  droplet  diameter 
increases,  the  injection  spray  mass  density  required  to  produce  the  same  amount  of  suppression 
also  increases.  Thus  in  a  co-flow  configuration,  smaller  diameter  mist  droplets  result  in  better  and 
more  efficient  suppression  for  the  same  injection  mass  density.  The  profile  of  peak  temperature  as  a 
function  of  injection  number  density  (Figure  23)  show  that  the  trends  are  reversed.  These  profiles 
show  that  larger  diameter  droplets  produce  more  effective  suppression  for  a  fixed  inlet  injection 
number  density.  This  is  because  when  the  injection  number  density  is  held  fixed,  larger  diameter 
droplets  result  in  larger  amount  of  mass  entrained  into  the  flame  as  compared  to  smaller  diameter 
droplets. 


Figure  24  shows  the  centerline  H2O  gas  density  for  various  spray  injection  densities  of  50fi 
droplets.  This  figure  also  shows  the  base  case  centerline  H2O  gas  density  profile.  Since  water  mist 
is  absent  in  the  base  case  calculation,  this  density  value  is  purely  due  to  convection  and  diffusion  of 
water-vapor  produced  during  chemical  reactions  alone.  The  difference  between  the  spray  injection 
density  profile  and  the  base  case  profile  gives  a  rough  measure  of  the  amount  of  H2O  gas  entrained 
into  the  flame.  As  expected,  it  is  found  that  as  the  spray  injection  density  increases  the  amount  of 
mist  entrained  into  the  flame  also  increases  for  a  given  initial  droplet  diameter.  Figure  25  shows 
the  centerline  temperature  profiles  for  various  spray  injection  densities.  The  droplet  diameter 
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and  droplet  velocity  are  the  same  as  for  the  case  shown  in  Figure  24.  As  the  droplet  density 
increases,  the  centerline  temperature  reduces  and  this  difference  in  temperature  represents  the  net 
suppression  due  to  a  combined  effect  of  thermal  cooling  and  oxygen  dilution.  It  should  be  noted 
that  temperature  reduction  at  the  centerline  is  negligibly  small  during  the  first  1.5  cm  above  the 
burner  surface,  but  then  increases  with  the  height.  Figure  26  shows  the  amount  of  mist  entrained 
as  a  function  of  water  mist  droplet  diameter.  The  figure  clearly  shows  that  at  small  heights  (10mm) 
the  amount  of  mist  entrained  reduces  as  the  droplet  diameter  increases.  Larger  diameter  droplets 
preserve  their  initial  momentum  for  a  longer  time  than  the  smaller  diameter  droplets.  As  a  result, 
the  larger  diameter  droplets  take  a  longer  time  to  be  entrained  into  the  flame.  Figure  27  shows  the 
amount  of  mist  entrained  at  the  centerline  as  a  function  of  spray  injection  flux  density  for  three 
initial  droplet  diameters  and  at  three  different  heights  above  the  burner  surface.  We  observe  that 
the  amount  of  entrainment  increases  linearly  with  spray  injection  flux  density  (mjfjo/^F)-  Droplet 
entrainment  increases  with  height  above  the  burner  surface  for  a  given  droplet  diameter.  Figure  27 
further  shows  the  impact  of  droplet  diameter  on  entrainment  rate  above  the  burner  surface.  It  is 
clearly  shown  that  the  smaller  diameter  droplets  entrain  more  readily  into  the  flame  whereas  larger 
diameter  droplets  entrain  at  a  smaller  rate,  for  a  given  spray  injection  density  and  at  a  given  height 
above  the  burner  surface.  This  result  is  strictly  valid  only  when  the  mist  is  introduced  parallel  to 
the  air-flow. 

The  net  suppression  due  to  water  mist  and  the  impact  of  droplet  diameter,  injection  velocity 
and  spray  density  have  been  summarized  in  Figure  28.  The  abcissa  of  this  figure  shows  the  net 
amount  of  water  mist  injected  per  unit  mass  flow  rate  of  fuel,  whereas  the  ordinate  is  a  ratio  of 
the  integrated  heat  release  in  the  computational  domain  with  water  mist  (AT^t)  to  the  integrated 
value  without  water  mist  (Affo)-  This  ratio  therefore  represents  a  net  integrated  measure  of 
the  suppression  of  a  flame.  Efficient  design  of  water  mist  systems  aims  at  obtaining  the  maximum 
amount  of  suppression  with  the  minimum  amount  of  water  mist  added  to  system,  that  is  to  minimize 
the  value  of  the  abcissa  and  the  ordinate.  Results  for  various  droplet  diameter  (50,  100  and  150//) 
with  various  injection  velocities  show  the  net  suppression  that  is  observed  for  the  various  cases. 
The  ratio  of  the  water  supply  rate  to  the  fuel  rate  required  for  extinction  of  the  diffusion  flame  is 
approximately  10  for  50//  droplets  and  increases  to  40  for  150//  droplets.  For  each  droplet  diameter, 
the  net  suppression  increases  with  spray  density.  Figure  28  also  shows  three  cases  in  which  the 
spray  injection  velocity  is  increased  from  25cm/sec  to  500c7n/sec.  The  droplet  diameter  for  these 
cases  is  150//.  As  the  spray  velocity  increases,  we  find  that  net  suppression  reduces  for  a  given 
mist  flux  density.  Overall  our  results  indicate  that  for  the  co-flow  configuration,  smaller  diameter 
droplets  produce  maximum  suppression  for  the  minimum  spray  mass  density. 

4.  CONCLUSIONS 

A  numerical  model  has  been  developed  to  study  the  combustion  of  methane-air  diffusion  flames 
and  their  inhibition  by  water  mist.  Both  the  gas  phase  and  the  droplet  phase  flows  are  described  by 
equations  in  the  Eulerian  form.  The  complete  set  of  unsteady  compressible  Navier  Stokes  equations 
are  solved  for  studying  methane  -air  diffusion  flames  stabilized  above  Wolfhard  Parker  burners.  A 
sectional  approach  has  been  adopted  in  which  the  entire  droplet  size  domain  is  divided  into  a 
few  sections  and  an  integral  quantity  representing  each  section  is  conserved.  The  model  has  been 
validated  by  comparing  temperature  profiles  with  experimental  measurements.  The  model  has  been 
used  to  understand  the  impact  of  droplet  diameter,  spray  velocity  and  injection  characteristics  on 
mist  entrainment  into  the  diffusion  flame  and  flame  suppression.  Numerical  simulations  have  been 
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performed  to  identify  the  relative  contribution  of  the  various  suppression  mechanisms.  Our  results 
indicate  that  smaller  droplets  exhibit  smaller  characteristic  time  for  decrease  in  relative  velocity, 
with  respect  to  the  gas  phase,  and  therefore  entrain  most  rapidly  into  the  diffusion  flame.  Hence, 
for  the  co-flow  injection  case,  smaller  diameter  droplets  produce  maximum  flame  suppression  for 
a  fixed  amount  of  injection  spray  density.  The  ratio  of  water  supply  rate  to  the  fuel  supply  rate 
required  for  extinguishment,  reduces  with  decreasing  droplet  size  for  laboratory  flames. 

Additional  studies  are  being  pursued  with  different  mist  injection  configurations  such  as  side 
injection,  top  injection  or  injection  at  an  angle  to  the  co-flow  to  further  optimize  the  various  water 
mist  injection  characteristics  for  flame  suppression.  The  net  suppression  (AFr)  for  the  various 
mist  injection  configurations  will  be  computed  and  presented  in  a  form  similar  to  the  one  shown 
in  Figure  28.  The  suppression  of  non-symmetric  flames  when  subjected  to  side  injection  of  fine 
water  mist  is  being  investigated  using  numerical  simulations  in  which  the  symmetry  condition  is 
no  longer  enforced.  These  studies  would  help  in  optimizing  spray  injection  characteristics  for  flame 
suppression.  Another  part  of  this  research  is  directed  toward  simulating  methanol  or  heptane 
liquid  pool  fires  and  simulating  the  changes  in  burning  rates  when  subjected  to  fine  water  mist 
droplets.  The  effect  of  droplet  diameter,  spray  injection  velocity  and  mist  density  on  methanol 
burning  rates  will  be  investigated.  The  effect  of  gas  velocity  on  droplet  entrainment  into  the 
diffusion  flame  needs  to  be  studied  further.  Instead  of  the  eulerian  approach  used  for  studying 
droplet  evaporation,  a  lagrangian  approach  for  tracking  the  individual  droplets  might  also  yield 
useful  information  regarding  the  movement  of  the  droplets  in  a  diffusion  flame  flow  field.  New 
improved  models  for  studying  the  radiation  attenuation  of  fine  water  mist  are  also  required. 
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Fig.  1.  Schematic  diagram  of  a  2-Dimensional  Wolfhard  -  Parker  diffusion  flame  burner,  showing 
burner  geometry.  The  burner  has  been  modified  to  introduce  mist  into  the  air  co-flow  through  mist 
generation  chambers  located  next  to  the  oxidizer  channel.  A  fraction  of  the  mist  generated  in  this 
chamber  is  entrained  into  the  air  stream  and  interacts  with  the  diffusion  flame. 
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Fig.  2.  Comparison  of  numerical  and  experimental  results  for  base  case  CH4  —  Air  diffusion 
flame.  Numerically  computed  temperature  profiles  are  compared  with  thermocouple  temperature 
measurements  at  a  height  of  5,  15,  30  and  55  mm  above  the  burner  surface. 
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Fig.  3.  Comparison  of  numerical  and  experimental  results  for  CH4  and  18.2%02,81.8%iV2  diffusion 
flame.  Numerically  computed  temperature  profiles  are  compared  with  thermocouple  temperature 
measurements  at  a  height  of  6,  40,  46  and  50  mm  above  the  burner  surface. 
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Fig.  4.  Comparison  of  numerical  and  experimental  results  for  CH4  and  16.5%025  83.5%iV2  diffusion 
flame.  Numerically  computed  temperature  profiles  are  compared  with  thermocouple  temperature 
measurements  at  a  height  of  6,  26,  46  and  60  mm  above  the  burner  surface. 
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Fig.  5.  Contours  of  temperature  and  heat  release  rate  over  a  CH4,  -  Air  diffusion  flame  burner.  The 
fuel  jet  velocity  is  2.81cm/ sec  and  air  jet  velocity  is  25.87cm/sec.  The  flame  height  is  approximately 
7cm  and  the  maximum  flame  temperature  is  1970A'.  The  heat  release  rate  profile  shows  the  region 
of  exothermic  chemical  activity  above  the  burner  surface. 
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Fig.  6.  Contours  of  methane  CHi  and  oxygen  O2  density  above  a  —  Air  diffusion  flame 
burner.  The  fuel  jet  velocity  is  2.%lcm/sec  and  air  jet  velocity  is  25.87cm/sec.  These  calculations 
are  performed  using  a  single  step  arrhenius  chemistry  to  study  methane-air  combustion. 
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Fig.  8.  Velocity  vectors  and  stream  line  pattern  during  methane-air  combustion  above  a  Wolfhard 
Parker  dilFusion  flame  burner.  The  velocity  vectors  have  been  color  coded  with  red  showing  the 
highest  gas  speed  and  blue  showing  the  lowest  gas  velocity.  The  streamline  pattern  (black)  clearly 
indicate  the  entrainment  of  oxidizer  along  with  the  suppressants  into  the  diffusion  flame.  The  fuel 
jet  velocity  is  2.81cm/sec  and  air  jet  velocity  is  25.87cm/sec. 
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Fig.  9.  Suppression  of  diffusion  flames  subjected  to  50/i  diameter  fine  water  mist  droplets  injected 
with  a  velocity  of  25cm/ sec  into  the  air  flow.  Figure  shows  reduction  in  flame  temperature  and  heat 
release  rate  profiles.  The  temperature  and  heat  release  rate  profiles  indicate  an  increase  in  flame 
height  and  flame  suppression  due  to  sub-extinguishing  application  rates  of  water  mist.  The  heat 
release  rate  profile  also  indicates  the  presence  of  an  endothermic  heat  release  profile  (-ve  values) 
showing  the  approximate  location  of  evaporation  of  the  droplets. 
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Fig.  11.  Suppression  of  dilfusion  flames  subjected  to  50/i  water  mist  droplets.  Figure  shows 
temperature  contours  for  various  spray  injection  densities  of  1000, 2000, 3000  and  4000drops/cm®. 
These  droplets  are  injected  along  with  the  air  with  an  initial  injection  velocity  of  25cm/sec.  Also 
shown  are  the  base  case  profiles  indicating  various  levels  of  suppression. 

The  flame  height  increases  from  7cm  for  the  base  case  to  11cm  for  injection  of  dOOOdrops/cm^  and 
the  flame  temperature  decreases  from  1970A'  to  1695/ir. 
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Fig.  12.  Sectional  droplet  density  contours  for  the  injection  of  50 fi  water  mist  droplets  along  with 
the  air  co-flow.  The  droplets  are  divided  into  5  sections  as  follows:  0  —  10/i,  10  —  20fi,  20  —  30/t, 
30  -  40^  and  40  -  50/i.  Droplet  section  densities  for  three  of  the  five  sections  have  been  shown 
to  illustrate  the  location  of  maximum  density  and  position  where  these  droplets  evaporate  to  form 
smaller  and  smaller  droplets  (sections). 
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Fig,  14.  Velocity  profiles  of  the  various  sections  have  been  shown  as  a  function  of  distance  above  the 
burner  surface  for  the  injection  of  50^  water  droplets.  The  velocity  component  shown  is  normal  to 
the  burner  surface.  The  droplets  are  divided  into  5  sections  as  follows:  0  30^,  30  —  60/x,  60  —  90^, 

90  -  120/x  and  120  -  150/i. 

Also  shown  are  the  gas  phase  profiles  indicating  the  effect  of  droplet  drag  and  buoyancy  forces  on 
various  droplet  sections.  This  figure  illustrates  the  different  characteristic  times  for  reduction  of 
relative  velocity  for  the  various  droplet  diameters. 
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Fig.  15.  Comparison  of  numerically  computed  and  experimentally  measured  temperature  profiles 
at  a  height  of  10mm,  30mm  and  45mm  above  the  burner  surface  with  injection  of  3%  mass  fraction 
steam  along  with  the  air  co-flow. 
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Fig.  16.  Comparison  of  numerically  computed  and  experimentally  measured  temperature  profiles 
at  a  height  of  10mm,  30mm  and  45mm  above  a  burner  surface  for  injection  of  3%  mass  fraction 
water  mist  droplet  of  50^  diameter. 
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b)  DiflFusion  Flame  with  Water  Mist 
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Fig.  17.  Experimentally  measured  temperature-time  history  at  two  different  points  (5mm, 30mm) 
and  (0mm,10ram)  for  the  a)  base  case  methane-air  diffusion  flame,  b)  diflfusion  flame  subjected  to 
3%  mass  fraction,  50p  water  mist  droplets  injected  along  with  the  air  co-flow.  Figure  illustrates 
that  the  diffusion  flame  subjected  to  water  mist  fluctuates  rapidly  and  is  shown  by  the  changes  in 
thermocouple  temperature  measurements. 
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Fig.  18.  Comparison  of  numerically  predicted  and  experimentally  measured  temperature  profiles  as 
a  function  of  distance  above  the  burner  surface  for  the  interaction  of  3%  mass  fraction,  50/x  water 
mist  droplets  with  the  diffusion  flame.  The  numerical  result  is  the  maximum  temperature  measured 
at  a  given  height  above  the  burner  surface.  The  experimentally  measured  maximum,  mean  and 
minimum  temperatures  are  the  maximum  values  at  a  given  height. 
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Fig.  19.  Temperature  contours  above  a  methane-air  diffusion  flame,  a)  Base  case  contours,  b) 
Temperature  contours  during  the  injection  of  50/x  water  droplets  with  a  spray  injection  density  of 
2000  drops/ c)  Same  as  case  b)  but  with  latent  heat  of  evaporation  set  to  zero.  These  contours 
illustrate  the  effect  of  the  various  suppression  mechanisms. 
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Fig.  20.  Temperature  profiles  at  various  heights  above  the  burner  surface  (for  the  three  cases 
discussed  in  the  previous  figure,  showing  the  relative  contribution  of  the  various  suppression 
mechanisms  such  as  oxygen  dilution  and  the  thermal  cooling.  Water  mist  droplets  are  50/l/  diameter. 
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Fig.  21.  Temperature  profiles  at  various  heights  above  the  burner  surface  (for  the  three  cases 
discussed  in  the  previous  figure,  showing  the  relative  contribution  of  the  various  suppression 
mechanisms  such  as  oxygen  dilution  and  the  thermal  cooling.  Water  mist  droplets  are  150^  in 
diameter. 
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Fig.  22.  Suppression /extinction  of  methane-air  diffusion  flames  during  injection  of  water  mist  along 
with  the  air  co-flow.  Figure  shows  the  reduction  in  maximum  flame  temperature  as  a  function  of 
spray  mass  density  for  various  droplet  diameters. 


Fig.  23.  Suppression/extinction  of  Methane- Air  diffusion  flames  during  injection  of  water  mist  along 
with  the  air  co-flow.  Figure  shows  the  reduction  in  maximum  flame  temperature  as  a  function  of 
droplet  number  density  for  various  droplet  diameters. 
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Fig.  24.  Center  line  H2O  gas  density  profile  as  a  function  of  distance  above  the  burner  surface  for 
various  spray  injection  densities  of  50 fi  water  droplets.  Also  shown  is  the  base  case  profile  indicating 
the  degree  of  entrainment  of  water  mist  into  the  flame. 
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Fig.  25.  Center  line  temperature  profile  as  a  function  of  distance  above  the  burner  surface  for  various 
spray  injection  densities  of  50//  water  droplets.  Also  shown  is  the  base  case  profile  indicating  the 
suppression  in  flame  temperature  due  to  injection  of  water  mist. 
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Fig.  26.  Effect  of  droplet  diameter  on  water  mist  entrainment. 
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Fig.  27.  Effect  of  droplet  diameter  and  injection  spray  density  on  water  mist  entrainment. 
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Fig.  28.  Integrated  effect  of  water  mist  on  diffusion  flames.  The  figure  shows  the  net  suppression 
effect  of  water  mist  and  its  dependence  on  injection  spray  density  and  velocity.  The  effect  of  droplet 
diameter  on  suppression  has  also  been  shown. 


